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~^ Abstract. We introduce an algorithm that iteratively produces a sequence of nat- 

I ural numbers ki and functions bi defined in the interval [l,+oo). The number /jj+i 

^ arises as the first point of discontinuity of bi above ki. We derive a set of properties 

_i of both sequences, suggesting that (1) the algorithm produces square- free numbers 

psj ki, (2) all the square-free numbers are generated as the output of the algorithm, and 

^_^ (3) the value of the Mobius function fx{ki) can be evaluated as 6i(/cj+i) — bi{ki). The 

^ logical equivalence of these properties is rigorously proved. The question remains 

^^ open if one of these properties can be derived from the definition of the algorithm. 

r^ Numerical evidence, limited to 5 x 10^, seems to support this conjecture. 

P Keywords: Mobius function, square-free numbers, zeta function, Riemann 
hypothesis. 
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in 

^- 1 Introduction 

o 

T-H A natural number n G N is called square-free, if the exponents arising in its prime factor- 

T^ ization 

"^ n = p^^ p^2 ' ' ' Pk 

are all equal to 1, i.e., ri = r2 = ■ ■ ■ = r^ = 1. For a natural number n with prime 



X 

^ factorization as above, the Mobius function jj, is defined as 



1 


ii n = 1; 


(_l)Ei.n, 


if n is square-free; 


0, 


otherwise. 



^{n) 



In other words, ^{n) is zero when n has a square factor, and otherwise gives the parity of 
the number of (distinct) prime factors of n. The Mobius function has important applica- 
tions in number theory, many of them concerning to the Riemann hypothesis about the 
zeros of the zeta function, P, [6], [7], [8j. 



An outstanding problem in algorithmic number theory is to compute fi{n) efficiently 
without first factoring n. By "efficiently" we mean a number of bit operations bounded by 
a polynomial in logn, the length of n in binary. As far as we know, the question remains 
open if the computation of fi{n) can be done in polynomial time, and in fact, nobody 
currently knows a way to compute it significantly faster than factoring n. 

In this paper, we present an algorithm that iteratively produces a sequence of numbers 
ki and the value of /i(fcj). In order to determine ^{ki), it is necessary to generate the whole 
sequence ki,k2, . . . ,ki. Our algorithm is based on a sequence of arithmetical functions bi, 
with the numbers ki arising as discontinuity points. These functions are closely related to 
the Nyman-Beurling approach to the Riemann hypothesis. The main result of the present 
paper is a set of properties of the sequences ki and bi, suggesting that 

• The numbers ki generated by the algorithm are square- free. 

• The set of all the square-free numbers can be generated by the algorithm. 

• The value of the Mobius function fi{ki) can be evaluated as bi{ki+i) — bi{ki). 

We are able to prove the logical equivalence of these properties. Unfortunately, using the 
definition of the algorithm, we cannot prove, neither disprove, if one (then all) of these 
conditions are satisfied. Numerical evidence, hmited to 5 x 10^, seems to support the 
conjectures quoted above. 

This paper is organized as follows. The section [2] introduces the framework of the 
Nyman-Beurling approach to Riemann hypothesis. The section |3] provides the basic def- 
initions for this paper, and the algorithm is defined in section |4| The main result of the 
present work is stated and proved in section [5} while further remarks are done in section 
|6} Finally, the implementation of the algorithm and numerical results are discussed in 
section [71 

2 Hilbert Space Approach to Riemann Hypothesis 

Denote by [a;J the integer part of x, i.e., the greatest integer less than, or equal to, x. 
Define the fractional part function by {x} = x — [x\. Given n G N and two families of 
parameters {ak}^^i C C and {Ok}^=i C (0, 1], we define a Beurling function as a function 
Fn (the sub-index n included in the notation for convenience) of the form 

Fn{x):=J2<^JH. (1) 

fc=l *- ^ 



For a Beurling function F„, an elementary computation shows that 

n 



{Fn{x) + l)x'~^dx 




c(^)E«'^^^ 



(2) 



A;=l 



for the complex variable s G C in the half-plane Re(s) > 0, i.e. for s with positive real 
part. Here, ( denotes the Riemann's zeta function given by 



OO -. 



fc=l 



The classic references are [8] and [6]. A derivation of the relation (pi) can be found, for 
instance, in [21 p. 253]. It is useful (but not always necessary from a theoretical point 
of view) to assume that the parameters defining the function F„ satisfy the additional 
condition 



7 _, O'k Sk 



0. 



(3) 



fc=i 



In this case, the first term at the right-hand side of ^ vanishes, simplifying the expression. 
The identity ^ is the starting point of the following theorem by Beurling 

Theorem 1 (Beurling). The zeta function C(s) has no zeros in the half-plane Re(s) > 
1/p if and only if the set of (Beurling) functions {fe{x) = {0/x}}o^g<^i is dense in 
LP{[0,l],dx). 

See [21 p. 252] for a proof of the Beurling theorem and further references. Note that for 
p = 2 this result provides an equivalent condition to the Riemann hypothesis (RH) for the 
zeta function. This is the Beurling, or Nyman- Beurling, approach to RH. 

The Beurling theorem above has an easy half part, whose proof can be sketched as 
follows. From relation (pi) and assuming ([3]), we have 



as)J2a,9l 



fc=i 



(F„(x) + l)x''~^dx 



^ \\Fn{x) + l\ 



„s-l| 



\X 



where the last relation follows using the Schwarz inequality in L'^{[0,l],dx). Therefore, 
if the first norm in the right-hand side above can be done arbitrarily small for a suitable 
choice of n, a^'s and O^s, then the function ({s) could not have zeros for Re(s) > 1/2. We 
will refer to the first condition above as the Beurling criterion (BC) for RH. Note also that 
in order to demonstrate the RH, it is sufficient to prove that the constant function equal 
to —1 can be arbitrarily approximated in the norm of the Hilbert space L^([0, 1], dx) by 



Beurling functions F„ of the form ([T|. It was proved in P] that BC remains equivalent to 
RH if the parameters Ok are restricted to be reciprocal of natural numbers, i.e. 9k = 1/bk, 
with bk e N. 

Several approximating functions to —1 of the form ([I]) were proposed in the literature. 
From the relations ^ and ([3]), we have that under the BC, the "partial sum" 

n 

Y^akOl (4) 

k=l 

is an approximation to the reciprocal of the zeta function l/(^(s), which is known to have 
an expression as a Dirichlet series 



C(») i:r * 



fc=l 



convergent for Re(s) > 1. Therefore, a (naive) first choice for an approximating function 
would be 



k=i I ^ ) 



(6) 



Note that this function does not matches the condition ^. We can handle this without 
subtlety, just by subtracting the difference, that is given by g{n), where 

N3k^t 

Therefore, a second choice would be 

B„(x) = X]^(*) i^\ - ng(n) |^} (8) 

g.w{f} -««(.-!) {i^}. (9) 



k=l 

Other variants were also proposed, as 



Vn{x):=±m[^]-9in)!^l^ 



(10) 



Unfortunately, the sequences (|6]), ^ and (10) are known to be not convergent to —1 
in L'^{[0, 1], dx), as proved in ^. A survey on the Nyman-Beurling reformulation of the 
Riemann hypothesis and later developments by Baez-Duarte can be found in |5j. 



3 Basic Definitions 

In order to motivate the definitions below, assume that a Beurhng function Fn as in ([I]) is 
constant between the reciprocal of the natural numbers. In other words, assume that such 
a function takes a constant value in each of the intervals (^rj, |], for all /c G N (but the 
constant value may differ from interval to interval). In this case, the integral in pi) can 
be expressed alternatively as an infinite series, involving the values of fn{k) ■■=Fn{l/k). 
Furthermore, the values of Fn{x) for all x G [0, 1] are completely determined by the values 
of fn{k) for all /c e N. We will call an arithmetical Beurling function a function of the 
form fn{k) = Fn{l/k), where F^ is a Beurling function. 

We introduce now, perhaps the simplest, non-trivial, example of arithmetical Beurling 
function satisfying the condition ([s]). For a,b eN define the function /3a,b as 

M^y- {-a} - 'a {^ ' (") 

As the functions Pa,b will be the basic blocks in our construction, we summarize some of 
its elementary properties in the following result. 

Lemma 1. Consider a,b eM., with < a < b. Then, 

a. < — > and \ -r\ are right- continuous, and linearly independent functions. 

b. (3a^b{x) = 0, when ^ x < a. 

c. Let k eN be such that {k — l)a < b ^ ka. Then, 

I ^ ^ j-j a ja^x < (j + l)a, for j = 1, . . . , (A; - 2); 

''''''^'''~\-(k-l) ii{k-l)a<:x<b. 

d. Assume a, 6 G N. Then, Pa,b{x) is constant when k ^ x < k + 1, for all /c G N. 



Proof: (a): The right-continuity is derived from of [x\. Now, ifci{|} -|- C2{|} = 
for all X, then for x = a we have = ci {^} + C2 {f}=C2f- Thus, C2 = and we have 
ci {^} = for all X, and taking now x = a/2 we get = ci |^} = ci/2 and ci = 0. 

(b): If ^ x < a < 6, then x/b < 1 and x/a < 1. Thus, {|} - ^ {|} = ^ - ^| = 0. 

(c): Assume j = 1, . . . ,{k — 2). Then, for ja ^ x < {j + l)a < b, we have x/b < 1 and 

J ^ x/a < (j + 1). Thus, {f} - Hf} = f - [fJ - !f = f - J - f = -J- Analo- 
gously, for {k — l)a ^ x < b < ka, we have x/b < 1 and (/c — 1) ^ x/a < k. Thus, 



S^X— b(x\_x_\x\_bx_x_^j^_'^\_x__^j^_'^\ 
\ a ) a \ b ) a La J a b a ^ ' a ^ ' ' 

(d): \i X < h then (d) is true by (b) and (c) already proven. Consider now b ^ x. If 
k < X < k + 1 then x ^ N and thus x/a ^ N and x/b ^ N. Therefore, there exists 
/cq and /q in N such that ak^ < x < a{ko + 1) and Mq < x < 6(/o + 1), and we have 
{a}~a{f} ~ a~^o~a(f~^o) = ~^o + ^^0; which is a coustaut independent of x. As 
{|} — ^{|}is right-continuous, by part (a), this is also true for A; ^ x < A; + 1. I 

4 The Algorithm 

We will define a sequence of numbers {fcijigN and functions {6i}jgN iteratively as follows. 
Start with the following definitions 



ki :=1 

k2 --=2 



''^<^>-^f}-|{l 



(12) 



Assuming now that ki and bi are already defined, for i ^ 2 define iteratively /cj+i and bi+i 
as follows. The number fcj+i is defined as /cj+i — ki + j, where j is the least integer such 
that bi{ki + j) 7^ bi{ki). Once determined the number /cj+i, the function bi+i is defined as 

W.)=*.W + (l + *.fe))({^}-^{^}). (13) 

Some elementary properties derived from these definitions are summarized in the following 
result. 

Lemma 2. For any i eN we have 

a. bi is a right- continuous function, which is constant between the natural numbers. 

b. bi+i{ki) = -1. 

c. Assume /cj+i ^ 2ki for i ^ 2. Then, bi{x) = —1 for all x G [1, ki). In particular, the 
sequence {6j}jeN converges point-wise to —1 in [l,+oo). 



Proof: (a): Observe that each bi is a (finite) hnear combination of f3p^q. Therefore, this 
resuh is a direct consequence of Lemma [l] 



(b): Is an immediate consequence of definition (13) 



(c): For induction on i. The case z = 2 is an immediate consequence of definition (12) 



Assume now that bj{x) = — 1 when x G [1, fcj) for all j ^ i. li x < ki < /cj+i, then from 
definition (13) we have 6j+i(x) = bi{x) which is equal to —1 by the inductive hypothesis. 



Now ii ki ^ X < fcj+i ^ 2ki, also from definition (|13|) we have 

■■ bi{x) - 1 - bi{ki) 



bi+i{x) =bi{x) + {l + bi{ki)) 



X 



1 - 



fcj-l-l X 



ki ki fcj+i 

because by the definition of fcj+i, the function bi{x) is constant for x G [ki, ki^i) 



I 



Remark: As bi is a constant function between the natural numbers, the number fcj+i is 
the first point of discontinuity of bi above ki. 



5 Main Result 

The next result is relevant in order to establish a relationship between the sequence 
{fcjjjgN, the values of /i(/ci), and the square-free numbers. 

Lemma 3. The following conditions are equivalent 
^{kj) _ 1 + bi{ki) 



-■E 



1 kj 



ki 



-, for i ^ 2. 



b.,w^x:My{^}-*.(i:^){|}./°"^2, 

fi(ki) l + bi{ki) l + bi-i{ki^i) 
c. — ^ — - = —^ — - — ; — ^ -, for 1^3. 



Furthermore, if the condition fcj+i < 2ki is valid for i ^ 2, then all conditions above are 
also equivalent to the following ones 

d. fi{ki+i) = bi{ki+i) - bi{ki), fori ^ 2. 
ki 



e. ^^l{kj 



L^JJ 



1, fori ^ 1. 



Proof: In order to prove the logical equivalence between all conditions in Lemma |3| we 
separately will prove, first of all, the equivalences (a)-v4>(b) and (a)<;=^(c). Then, after the 
introduction of the additional condition fcj+i < 2fcj, we will prove the equivalence between 
(d)<^(c) and (e)<^(b). 



(a)^(b): For induction on i. For i = 2, from definition (12) we have 
Assuming now that condition (b) follows for all j such that j ^ i, we have 



b,+i{x) = bi{x) + {1 + bi{h)) 



rCi I rCi I rC' 



H+l 



X \ I ^ — ^ jJjirCjj \l^| /-, T/7\\/l'^l "^i+l I ^ 



X \ 1 I \ — ^ /^V 7/ I I '^ 1 / 7 \ — ^ l~^\j) I / I "^ 1 '^i-\-l I ^ 



|.(.){a-J|f^){^} 



hjj^ J '^i I '^' 



^i+l 



i-1 ^ ^ /i-1 



j=l y J ^ \j=\ 






E.(y^f|-- E^)i^[;a5) 



and this proves condition (b) for i + 1. Here, in the first equality we have used the induc- 
tive hypothesis and in the second one we have used condition (a). 



(b)^(a): From Lemma p^(b), by condition (b) we have 
= l + bi+,iki) 



j=i ^ 1 -> \j=i 1 

i + EMyJI^ 



Ki 



i+1 



+ /i(fci) ]y}~ ^i+i 



fi{kj] 



E 



■^ kj I k,i+i 



(16) 



.i=i 



and this proves condition (a). Here we have used that \ ^^ \ = ^-^, (because -^-^ < 1), 

that ■) ^ f = 0, and also the relation bi{ki) = J2]=i f^i^j) i f^ f ' '^hich is an easy conse- 
quence of condition (b). 

(a)^(c): Using condition (a) for i and i — 1 we have 

fi{kj) _ 1 + bi{ki) 



E 



rC'i 

i-l 

^ k. 

j=l J 



K4 



(17) 



fJ'jkj) _ 1 + bi-i{ki_i) 
ki-i 



and subtracting (18) from (17) we get condition (c). 



(c)^(a): Denoting a{i) 






from (c) we have 



E 



H{kj) /i(A;i) iJ,{k2) 



j=i 



rv-i 



ki 



+ 



ko 



' 1 

+ ^a(j)-«(j-l) = 1 - ;^ + a(^)-a(2) 



i=3 



. 1 /. 1 + ^2(^2) 



1 - - + a(z) 



1 
2 



a[i] 



1 + &i(fc^ 

Ki 



(19) 



Here we have used that 62(^2) = 0, by definition (12). 

Therefore, we have proved (b)<(=^(a)<^(c). Assume now condition /cj+i < 2ki, for i ^ 2. 
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(d)<^(c): From definition (13) we have 

b.iki) = b,.i{ki) + (1 + bi_^{ki_i)) 
= 6,_i(/ci) + (l + 6i_i(A;,_i)) 
= 6,_i(/ci) + (l + 6i_i(fci_i)) 
= bi.iih) + {1 + bi^iih^i)) 



rh'i I rh'i 






ki^i 
ki-i 



ki-1 



(20) 



- 1 



b^-iih) + (1 + 6i-i(A;,_i)) ^ - (1 + 6i-i(A;._i)) . 



ki-i 



Observe that the additional condition implies 1 < -^-^ < 2, and therefore 



From (20) follows 



Ki- 



1. 



l + bi{ki) l + bi-i{ki_i) bi^i{ki) - bi-i{ki^i) 



Ki 



ki-1 



Ki 



(21) 



for i ^ 3. The equivalence between (d) and (c) is a direct consequence of (21) above. 
(b)^(e): From Lemma |2](b), by condition (b) we have 



= 1 +*.(*.-)= 1 + EMy {*|i} -«=. (i: ^j {*if} 



' i-l 



fi{kj 



i-l 



1+ Et=^ *--E''ft) 






'i-l 



+ «=. E 



fi{kj 



o=i 



k 


-1 




[k, J 




\ 


/ci-1 


/ 


[ h J 



i-l 



-*. E 



vj=l 



/i(fci) \ fc^-1 



i-l 






and this proves condition (e). Here, we have used that -r-^ < 1 and therefore 



ki 


-1 




[k, J 


; 


re 


k 


I 



(22) 



0. 



(e)^(b): Assume condition (e) valid for all i G N. We will prove (b) by induction on i. 



Condition (b) for i = 2 follows from definition (12) as done in the proof (a)^(b) above 

^i fi{l 



Assume now condition (b) valid for all n ^ i. Using the notation g{i) = ^^^^ h^^ ^^ 



10 



have 



j-i 






Ki 



i+1 



i-1 



*■ E 



fi{kj) \ r ki+i 



i-1 






fcj+1 



kig{i~ 1) ^±1 + /^.^(^_i) 



fcj+1 






Jy 



"-j+l 



+ kig{i- 1) 



(23) 



— I 1 — /i(A;. 



j+ly 



k 



i+1 



k 



i+1 



- fi{ki 



ki+i 



+ kig{i-l) 



Here we have used 
ogously we have 



- (1 - /i(/ci+i) - ^^{ki)) + kigii - 1) 
kig{i - 1) + ^i{ki+i) + /i(A;i) - 1. 



1 (a consequence of the additional condition) and (e). Anal- 



-1 



hi{ki) = y^/i(fc 



i=i "^^^ 



i-1 



MEt^ 



o=i 



E^(^j' 



i=i 



/Cj 



i-l 



kig{i- 1) -E^(^i) 



l'^jj 



(24) 



A;j5f(z - 1) - (1 - li{ki) 

kig{i - 1) - (1 - fi{ki)) 
kig{i- 1) + n{ki) - 1. 






Now, subtracting p4| ) from (23) we get condition (d). But we have already proved that 
(d)^(c)^(a)^(b). 



I 



6 Discussion 

The definitions in section |4] provide an algorithm to produce iteratively a sequence of 
numbers \ki\i<zfi and functions {6j}jgN- The Lemma |3] in section |5] states a circle of logi- 
cally equivalent properties of both sequences. 
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The condition (b) of Lemma p^ suggests that the functions 6j are the arithmetical coun- 
terpart of the approximating functions Bn in the relation ([9]). However, note that our 
definition of bi in section |4] is quite different of Bn- 



The condition (d) of Lemma p^ is related with the conjecture that the algorithm pro- 
duces square-free numbers fcj, and also provides the value of the Mobius function fi{ki). 
Note that the value of bi{ki^i) — bi{ki) is never zero, by the definition of the algorithm in 
section HI 



The condition (e) of Lemma M is related with the conjecture that the algorithm pro- 
duces all the square-free numbers. Indeed, if condition (e) were true, this would be a 
corollary of a well known result; see [Up. 66]. 



The condition fcj+i < 2ki in Lemma pi sufficient for (d) and (e), seems to be also 
necessary, as the following heuristic argument suggests. It is known that the square-free 
numbers are distributed in N with density 6/7r^; see [3 Thm. 333, p. 269]. Therefore, we 
can estimate the average distance between two consecutive square- free numbers as 7r^/6. 
Consequently, fcj+i ^ ki + 7r^/6, or ki+i/ki ;^ 1 -|- n'^/Gki. The last expression is less than 
2 for ki > 7r^/6 ~ 1.64. Thus, condition ki+i < 2ki for i ^ 2 seems to be reasonable also. 

Unfortunately, using the definition of the algorithm given in section |4], we cannot 
prove, neither disprove, if one of the conditions in Lemma |3] are satisfied. 

Note also that the algorithm cannot compute isolated values of fi{ki). In order to 
determine n{ki), it is necessary to generate the whole sequence ki,k2, ■ ■ ■ ,ki. 

7 Numerical Results 

The algorithm defined in section |4] was implemented using the Java programming lan- 



guage. The source code can be downloaded from http:// 143. 107.59. 106:96207cainille/ 
beurling.tar.bz2, This archive provides, in fact, two slightly different implementations 
of the algorithm. 

7.1 The class Beurling 

The methods in this class compute the sequences ki and bi, storing the values in an array 
with fixed size. The main method takes a natural number n as input. Its output is a file 
containing 3-uplas {ki, fj,{ki),ti), for i from 3 to n. Here, tj is the running time, in seconds, 
between the computation of /x(/ci_i) and fi{ki). The values of the Mobius function are 
calculated using the identity in Lemma [31(d). The main method in this class also verifies: 
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• If each one of the numbers ki is square- free. 

• If there exist eventually square-free numbers between ki and /c^+i that are not gen- 
erated by the algorithm. 

• The condition fcj+i < 2ki, equivalent to the gap fcj+i — ki < ki. 
These additional verifications are not included in the running time tj. 

Running this class with n = 1 x 10^, the size of the output file was about 10 MB. To 
estimate the running time, we generate a graphic with the pairs {ki,ti); see figure [1} 




266000 






/ 



























Figure 1: Pairs (ki,ti), where t, is the 
running time (seconds) between the com- 
putation of fi{ki_i) and fi{ki), using the 
class Beurling. 



Figure 2: Pairs (ki,Ti), where Tj is the 
total running time (seconds) for the com- 
putation of fi{ki), using the Java class 
Beurling. 



This graphic seems to suggest that square-free numbers can be divided in classes, and 
for each of these classes, the running time starting from the previous iteration t, is con- 
stant. We verify that all the square-free numbers in the range analyzed are generated 
by the program, and satisfy the condition fcj+i < 2ki, for z ^ 2. To estimate the total 
running time Tj to compute fi{ki), we use the formula Tj = ^^^i^fc- The graphic with 
the pairs (ki,Ti) is shown in figure 12} 



The program produced also four square-full, i.e. non square-free, numbers, shown 
in the table [Ij All these square-full numbers ki satisfy bi{ki+i) — bi(ki) = iJ,{ki) = 0. 
However, as already stated, by the definition of the algorithm in section |4} it must be 
bi{ki^i) — bi{ki) 7^ for any generated number. Therefore, we strongly believe that the 
square-full numbers are generated purely by rounding errors. 
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Number k 


^^{k) 


Divisible by 


440375 





52 


551208 





2^ 


799460 





2^ 


979275 





52 



Table 1: The four square-full numbers produced by the class Beurling, with ra = 1 x 10^ 



7.2 The class BeurlingArrayList 

This class works as the previous one, but the computed values are stored in an array with 
dynamic size. We use this class to process n = 5 x 10^ numbers, using a relatively modest 
desktop machine. The graphics with the pairs (A;j,tj) and {ki,Ti) are shown in figures [3 
and|4} respectively. 
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l.SetBB 2e+B6 2,geta6 3etaB 3^Se+B6 -le+BB I.SetaB Se+BB 



le+B6 Ii5e+B6 2e+B6 2,Seta6 3etaB 3,ge+BB 46+06 -l^Se+OB Se+BB 



Figure 3: Pairs (fci,tj), where tj is the 
running time (seconds) between the com- 
putation of /i(fci_i) and /i(fcj), using the 
class BeurlingArrayList. 



Figure 4: Pairs (fcj,Tj), where Tj is 
the total running time (seconds) for the 
computation of ^{ki), using the class 
BeurlingArrayList. 



As in the case of figure [TJ the square-free numbers ki seems to be divided in classes. 
However, in this case, the running time between two successive iterations tj grows linearly 
with ki- All the square- free numbers in the range analyzed are generated by the program, 
and satisfy the condition fcj+i < 2ki, for i ^ 2. In the range analized, twenty square-full 
numbers were also generated by this class, probably by rounding errors, as explained 
above. 
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8 Concluding Remarks 

In section |4] we define an algorithm that iteratively produces a sequence of numbers ki 
and functions 6j. The lemma [3] states a set of properties of these sequences suggesting 
that 

• The numbers ki generated by the algorithm are square- free. 

• The set of all the square-free numbers can be generated by the algorithm. 

• The value of the Mobius function /i(/ci) can be evaluated as fi{ki) = bi{ki^i) — hi{ki). 

In section [5] we prove the logical equivalence of these properties. Unfortunately, using the 
definition of the algorithm, we cannot prove, neither disprove, if one of these conditions 
are satisfied. Note also that in order to determine /i(fci), it is necessary to generate the 
whole sequence ki,k2, . . . ,ki. 

Numerical evidence seems to support the conjectures quoted above. However, this 
evidence is limited, and certainly not conclusive. 
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